Gross-Pitaevskii Equation 
for a System of Randomly Interacting Cold Bosons 
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Random interaction models have been successful in describing the amorphous properties of solids 
such as spin-glasses and structural glasses. This modelling approach is applied to a system of 
zero-spin cold bosons moving in an amorphous environment. The bosons are given to interact 
according to a non-random hard-core interaction. Additionally the bosons are subjected to a random 
interaction potential similar to that used for glasses. The approach is to apply a combination of 
replica methods and field theoretic techniques developed for superfluid Bose systems. This leads to 
an equation for the low-temperature Bose-Einstein condensate which is derived in the Hartree-Fock 
approximation. The equation is similar to the Gross-Pitaevskii equation, but the hard-core coupling 
constant is renormalised by the presence of the random interactions in a position dependent way. 
The amorphous contribution contains the replica diagonal and off-diagonal Green's functions, for 
which the Dyson equations are formulated. 



I. INTRODUCTION 

Random interaction potentials are used in models of 
amorphous materials like spin glasses ll| .quantum spin 
glasses [1, H, 0] and quantum glasses [3, 0| . Such mod- 
els are described by a system of interacting distinguish- 
able degrees of freedom. The present contribution is con- 
cerned with effects such interactions would have for a sys- 
tem of indistinguishable degrees of freedom. A system of 
cold bosons moving in an amorphous background with 
molecular-scale heterogeneities could be a candidate sys- 
tem described by such a theory if we think of the amor- 
phous background modulating the interactions. 

It is interesting to relate and contrast the present 
project with theories of Bose-glasses 0, H 111 describ- 
ing superfluid "'He confined to nano-porous glasses. The 
characteristic nano-scale of the disorder in these me- 
dia (the size of the pores) is considerably larger than 
the size of the bosons, and is typically described in 
terms of random external potentials. By contrast, in 
the present model the amorphous environment is repre- 
sented by random-interactions, and would more realisti- 
cally require molecular-scale heterogeneity. As a specific 
realization one might think of charged bosons moving 
in an environment with heterogeneous dielectric proper- 
ties, though it is an open question whether such system 
would be within reach of current fabrication techniques. 
For the time being our main focus will be on developing 
the theoretical description of such systems and to begin 
elucidating the effects random interactions could have. 

We shall formulate the theory in terms of coherent 
fields and a low-temperature coherent wave-function for 
which field-theoretic techniques are available, see e.g. 
[Tol . [Til . [T3 |. In particular this involves constructing 
a 'Gross-Pitaevskii like' equation for the Bose-Einstein 
condensate [l^ [l3| • The main aim here is to derive this 
equation and to establish the characteristic manner in 
which random interactions lead to a renormalisation of 
the bare (hard-core) coupling between the bosons. 

This paper is organised as follows. In section [ll] the 



random interaction Hamiltonian is introduced. Section 
mil gives the partition function in terms of a path integral 
in a coherent field representation and defines the disorder 
average. Then in section llVl the low-temperature conden- 
sate wave-function is introduced for which an equation is 
derived in section |V1 Its replica symmetric form is given 
in section IVTl Then the Dyson equations for the normal 
and anomalous Green's functions are given in section lVIII 
Finally, some conclusions are drawn in section rVIIII 



II. THE MODEL 

Starting point is the field-Hamiltonian of a three- 
dimensional system of interacting bosons with zero spin. 
Such requires the use of creation- and annihilation oper- 
ators [x) and V'l^c), formally creating and annihilating 
a particle at position x. The operators act on number 
states I {fix}) in a Fock space with variable number N 
of particles. The relations for these operators and the 
second quantization procedure to construct a Hamilto- 
nian are well known, see e.g. [Tsl . [l^ . In particular 
this assumes an interaction potential function V{x, x') 
for the interaction of a particle at position x with a par- 
ticle at position x' , where the vectors x and x' point 
to any position in the volume V. In what follows the 
three-dimensional integrals J dx are taken over the en- 
tire volume. 

We choose a random interacting potential of a similar 
type and interpretation as the ones used in mean-field 
glass and spin-glass theory. For the field theory pro- 
posed here this means that the amorphous background is 
represented by quenched random couplings J{x, x') de- 
scribing the interactions of particles at positions x and 
x'. We will take the J{x, x') to be zero mean Gaussians 
with J'^{x,x') = (j'^{x,x'). One would expect (t{x,x') 
to be translationally invariant, a{x,x') — a{\x — x'\), 
and decreasing sufficiently rapidly with distance. Alter- 
natively a mean-field model is constructed by choosing 
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(j'^{x,x') = /N , independently of distance \x — x'\. 

In addition to the random-interaction the bosons are 
given to interact according a hard-core interaction with 
(non-random) parameter Vo- Such approximation is jus- 



tified by the fact that at low temperatures the thermal 
de-Broglie wavelength of the bosons is much larger than 
the range of their interaction. In this setting we propose 
the following (grand canonical) field-Hamiltonian 



Hj[i,\i^] = j dx ^\x)(^- 

'^'^ J '^^ J 
+ — J dx J dx' 

where m is the mass of the bosons and fj, the chemical 
potential. The J-suffix in Hj emphasises the fact that 
this Hamiltonian represents a quenched amorphous back- 
ground characterised by the J{x,x'). The factor N is 
required to assure that the Hamiltonian is of 0{N). The 
grand partition function for the configuration J{x, x') is 

Zj = TieM-PHj[^\i;]). (2) 

In order to study the equilibrium properties of the 
model we need to evaluate the disorder averaged free en- 
ergy F. The 'replica trick' [T'| allows us to compute it 
as 



-/3F = logZj= hm-logZ^, (3) 

rwO n 

where n represents the number of replicas of the system, 
for which formally the limit zero is taken. Z" is the n 
times replicated partition function. The overline denotes 
the integrations over the Gaussians J{x^ x'), defining the 
average over all realizations of the random interaction. 




d{x - x') iPHx) iP^{x') ip{x') iPix) 
J{x,x')^Jj''{x)tp\x')i:{x') tPix), (1) 

I 

III. COHERENT FIELDS 



To obtain a path integral representation for the parti- 
tion function in the Matsubara formalism, we use a rep- 
resentation in terms of coherent states and (complex 
valued) coherent fields (j){x) [l^. In this representation 
the replicated partition function for the fixed disorder 
configuration J{x,x') is 



Zj = / n ^KT^'I'- exp( - - M])- (4) 



From here on we shall work in the scaled setting 
h — m = 1. The Euclidean action in the coherent field 
representation reads 



+ Y^J dT J (j)l{x,T)(j)l{x,T)<j)a{x,T)(f>a{x,T) 

+ dr J dx J dx' J{x,x') (j)l{x , t) (f>l{x' , t) (f>a{x' , t) (f>a{x , t) 



(5) 



r 



Next is to average ^ over the parameter J{x, x') by 
means of the following product of Gaussian integrations 



dJ{x, x') 



The notation Y[{x x') means that the integral is over pairs 
of points, assuming that J{x, x') is symmetric. The dis- 



exp 



(5) V2-^^(^>^') 

, J^x^x') 



dx / dx' 



2a^{x,x') 



(6) 
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order averaged replicated partition function becomes 

= y n T^'t^l'^^a exp( - 0J] ) . (7) 



The final Euclidean action S (for the disorder averaged 
theory) becomes 







dr / dx (I)1{x,t) (t>*Jx,T) (t)a{x,T) (t)a{x,T) 



4 ^ 



/3 /-/S 

dr / dr' 

Jo 



da; / dx' a'^{x, x') (t)l{x, t) (t)l{x', t) (t)a{x' , t) (t)a{x, t) 

X 0b (a;, r') (j3l{x' , t') <j)b{x', t') (t>b{x, t') 



(8) 



r 



In the coherent field representation the particle number 
constraint reads 

dx (0:(a;,r)0a(a;,r)} = N, (9) 
where (..) denotes the average w.r.t. the action S. 

IV. LOW-TEMPERATURE CONDENSATE 

Our focus is on the ultra-low temperature behaviour 
of the system, in particular the effect of Bose-Einstein 
condensation representing a dominant ground-state oc- 
cupation. This is taken into account by performing the 
shift 



(j)a{x,T) = «'a(a;) + (p'a{x,T), 



(10) 



in the functional integral ([7]) with the assumption that 
5'a(a;)| I d ) ' (x t)\ for typical configurations (j)'^{x,T) 
(see e.g. |10l.Tlll The shift 'ilaix) defines a non- 

fluctuating field and (j)'^{x,T) is the fluctuating field to 
be functionally integrated over. The partition function 
d?]) becomes 

/n 
1[V^':V4>: exp(-5*[{C,0l}]), (11) 
a=l 



where 5* denotes the shifted action 



(12) 



The form of 5* is complicated because a + <f>a appears 
to the powers four and eight in the interaction terms of 
([8|) . An approximated form of (I12p is treated in the next 
section. Averaging (fTO|) gives 



We require "ifaix) = (0a (a;, r))^* i.e. {(j)'^{x,T)) 3^ = 0. 
This is achieved by the condition that the linear terms in 
the shifted action vanish. Such condition leads to an 
equation for "ifaix) similar to the Gross-Pitacvskii equa- 
tion [1^ 13' derived in the next section. The particle 
number constraint ^ in the shifted theory becomes 



N = J dx |*a(a;)|' 
dx |*a(a;)|2. 



dx {<j)'*{x,T)<j)'^{x,T)) 



(14) 



V. HARTREE-FOCK APPROXIMATION 

Consider the path integral pT|) . The shifted action 5* 
may be approximated by a sum of an action containing 
only 4", an action term linear in 0' and an action term 
quadratic in cp' 



+ [{C0l;vl/:,vl;4] 

+ S^S^l[W:,^'a,K,^a}]. (15) 

The linear and quadratic term are constructed in the 
Hartree-Fock approximation. The quadratic part will be 
treated in section IVIII below. The linear action is 
required to vanish as mentioned in the previous section. 
It can be written in the general form 



c(HF) _ 

'^lin 



E 



dT j dx [<t>':{x,T) Va{x)<i>a{x) 
+ cp',(x,T)Va{x)n{x)]. (16) 



(0a(a;,r)}s* = ^a{x) + {(I)'^{x,t))s^. (13) 



The operator Va{x) is constructed from the action ([S]) 
with application of the shift (fTU)) . In the Hartree-Fock 
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approximation the result is 



b -^0 

_2 „/\iTr /~,/\iTr* /~,^lTr* /'~/\,T. /' ~, 



<a2(a;, x')^a{x')n{x)mx')^b{x), (17) 
where Gab{x'; t ~ t') denotes the Green's function 

Gab{x; T, t') = ( <P':{X, t) (t>'i,ix, t') ). (18) 

In an equihbriuni theory as presented here, the Green's 
function depends only on the time difference i.e. 
Gab{x; t,t') = Gab{x; t — t'). In pT]) we have exploited 
the fact that the spatially Ojff-diagonal Green's functions 
Gab{x, x'; T — t') vanish for all a, b [T^j. We have taken 
only the lowest order Hartree-Fock approximation. In the 
second term of we integrated by parts twice (w.r.t. 
V^) using periodic boundary conditions. Finally, in the 
amorphous part of ()17|) we omitted a term containing 
only the functions '^a{x),'i/a{x'),'ifb{x),'^b{x') and no 
Green's functions. This is justified by the fact that such 
term is 0{n) which vanishes in the replica limit n —> 0. 

To make sure that (0^(a;,r)} = in (fT3|) , the coeffi- 
cients of the linear term (|16p should vanish i.e. 

(19a) 
(19b) 



Va{x)^a{x) =0, 
Va{x)Kix) =0. 



The expression Vaix) defines a real operator since the 
quartic product of \l/-functions is real valued, so that we 
only need to consider (|19ap . We further assume that the 
Green's function Gab{x; r — r') is also real valued for all 
a, b. For reference we list here the complete equation 
(ITOa)) for ^'a(a;) 







b -^0 

xa^{x,x')^a{x')^l{x)^l{x')-i'b{x) 

^a{x). 



(20) 



This equation defines a 'Gross-Pitaevskii like' equation, 
with the contribution of the term in the second and third 
line representing the effect of the random interaction po- 
tential. In a replica symmetric approximation the struc- 
ture of this equation will become more transparent, which 
is subject of section IVll below. 

We also note that in deriving pT|) . use was made of 
the relation 



dx ^l{x) 



b',ix,T)+ J dx ^aix)<P':{ 



x,t)=0. (21) 



The physical reason behind this condition is that 0^(a;, t) 
should contain all the configurations that are orthogonal 

to ^aix) 0. 



VI. REPLICA SYMMETRY 

Next is to write down the replica symmetric (RS) equa- 
tion for the order parameter \E'(a;). At this early stage of 
developing the theory we do not consider replica symme- 
try breaking. The RS approximation is defined as 

^^J^{x) = ^{x), V(a) (22a) 
Gabix; t-t') = Gd{x; r - r'), V(a = b) (22b) 
Gab{x; t-t')^ G{x), V(a / b) (22c) 

Here we used the Ansatz that all replica off-diagonal 
Green's functions are time- independent. The argument 
to support this is that the replicas are independent and 
time-translation invariant, and that the origin of time 
could be chosen independently for each replica. In the 
RS approximation the equation (^0]) for \E'(£c) becomes 



= 



-\^{x)\^ dr' J dx' <j^{x,x') \^{x')\^ 

X ( Gd{x', t-t')- G{x') ) ] 'fix). (23) 



Gross-Pitaevskii equation and effective 
interactions 



Equation (j23p may be written in the familiar form of 
the Gross-Pitaevskii equation [l^, [l3| 







1. 



^i + Kff(/3,a;)|*(a;)|2 ^(^^^ (34) 



where now the coupling constant VcsiP, x) is an effective 
coupling constant defined as 



with 



Gdix') = - 

P Jo 



f3{Gd{x')-G{x')), (25) 



dr' Gd{x',T-T'), (26) 



This is one of our main results. Provided only that there 
is a non-zero condensate wave function ^(a;) 7^ the 
main effect of random interactions is a renormalisation of 
the bare (hard-core) interaction in a position dependent 
way. In the mean-field limit where a^{x,x') — cr^/iV, 
the renormalisation becomes position independent. 

Although typical values of the Green's functions are 
'small' - they are defined from small fluctuations (t>'{x, t) 
above the ground-state - we emphasize the fact that the 
function j3 {Gd{x') — G{x')) in the disorder contribution 
of ([25]) is not necessarily small, which is due to the ap- 
pearance of the factor (3. In particular at very low tem- 
peratures there could be a significant effect from this 
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factor. It is in this low-temperature region where also 
appearing in the amorphous term in (PS)) . is 
expected to become non-zero. 

In section IVlIl we shall formulate the Dyson equations 
for the Green's functions Gd and G. The Green's func- 
tions contain the wave- functions '^{x) as we shall see. 
A closed set of equations is obtained by adding to the 
Gross-Pitaevskii equation ([M]) and the Dyson equations 
the RS constraint equation (|14p reading 



TV 



j dx |*(a;)p + J dx Gd{x,0) w J dx |*(a;)p, 

(27) 

defining the normalisation condition for the condensate 
wave-function ^'(a;). 



and involve anomalous Green's functions. Anomalous 
Green's functions are a characteristic feature of theories 
involving expansions of the kind (fTO]) . see e.g. [13, [lH- 
We discuss here the general form of these equations. We 
shall work again in the Hartree-Fock approximation. As 
mentioned above, the normal Green's functions depend 
only on the time difference t — r' and are spatially diago- 
nal. This applies also to the anomalous Green's functions 
since they are determined from the normal Green's func- 
tions by the Dyson equations. 

In general the quadratic part of the shifted action (fTS]) 
may be represented in the terms of a Green's function 
matrix Qab{x] t — r') as 



VII. NORMAL AND ANOMALOUS GREEN'S 
FUNCTIONS 

The replica symmetric Dyson equations for the normal 
Green's functions Gd{x] t~t') and G{x) are complicated 



c(HF) r ,„ 



r 



The Green's function matrix reads 

(Gab{x;T ~ t') Aab{x;T - t')'' 
I '(29) 
Aab[x;T-T') Gab[x\T - t') 

where Gab are the normal and Aab and Aab the anoma- 
lous Green's functions. The latter are defined as 

Aab{x, T - r') ={ c^'^ix, t) <t,'^(x, t') ), (30a) 
Aab{x,r-T') ={c^':{x,t)c^'C[x,t')). (30b) 

It should be noted that the normal Green's functions are 
real valued but the anomalous Green's functions are not. 
The Green's functions satisfy the matrix valued Dyson 
equations [13,1111 

g-J{x-T-T')^g-l,{x-T^T')^i:ab{x-,T~T'), (31) 

with inverse propagator matrix 



Q^,!ab{^-^r^r')^G,^^,[x;r~r') (^J , 



(32) 



and inverse propagator 

d 1 



Golbi^-^ ^ - = - 2 V2 - 5ab 5{r - r'). (33) 



The self-energy matrix is expressed as 

'Y.<=^,{x:r-r') Y.t,{x:T - r'f 



T.ab{x;T - t') 



j:Ux-r-r') SG(a;;r-r'), 



.(34) 



Here the superscript G labels the normal self-energy and 
the superscripts A and A the anomalous self-energies. 
They are constructed from ([28]) . ((3T|) . ([34)1 and the action 
([S]) with application of the shift (fTU|) . In the Hartree-Fock 
approximation the result is 

-Y.^,{x:t-t') = 2VoK{x)^a{x)5ab5{T^T') 

-2 j dx' Gab{x'\T -T')a^{x,x') 

^■^a{x')-^a{x)mx)%{x') 

- j dx' Gab{x;T -T')a'^{x,x') 
x^l{x')^a{x')n{x')^,{x'), (35) 



^tf,{x-T-T') = Vo-^a[x)-^a{x) 5ab5{T-T') 

- j dx' Gab{x';T-T')a''(x,x') 

X ^a{x')^a{x)n{x')^b{x), (36) 
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-J dx' Gab{x'-T-T')a\x,x') 

x^i{x)-^a{x')n{^)n{x'). (37) 

Observe that the normal self-energies are real valued 
and the anomalous self-energies complex valued in gen- 
eral. The RS representations of the self-energies follow 
straightforwardly from (j22p . In a first approximation, 
simplifying the formulation considerably, we may ignore 
the amorphous parts in the self-energies since the Green's 
functions are 'small'. Finally, in the RS approxirnation 
the general form of the Dyson equations (pij) is Ij (in 
the replica limit n — > 0) 

{Qd - Q)-^ -Q{Qd- Q)-^ = Q^.li - Srf, (38a) 
-a(erf-e)-' = -S, (38b) 

where (|38ap defines the replica diagonal and (j38bp the 
replica off-diagonal part. These are matrix equations and 
constitute the six equations solving for the two normal 
Green's functions Gd, G and the four anomalous Green's 
functions Ad, A and Ad, A, each expressed in terms of the 
wave- function "^{x). 

As mentioned before, the precise solutions for ^'(a;) 
and the Green's functions are currently not yet known. 
When these become available one may compute the dis- 
persion relations for collective excitations from the poles 
of the ve^Yica,- diagonal Green's function. This amounts 
to evaluating in momentum space det {k; iuji) — 0, 
where uii are the Matsubara frequencies. 

VIII. CONCLUSIONS 

A random interaction model for a system of 
cold bosons is constructed, applying existing replica- 



techniques used in theories of spin-glasses and structural 
glasses. The random interaction was assumed to repre- 
sent an amorphous environment for the bosons. In ad- 
dition the bosons were given to interact according to a 
hard-core interaction. A field-Hamiltonian was proposed 
and a disorder averaged partition function in the coher- 
ent field representation was presented. A condition was 
derived for which a path integral shift was equal to the ex- 
pectation value of the field. This shift defined the wave- 
function for the Bose-Einstein condensate. The condition 
led to a replica symmetric equation for the wave- function 
which was similar to the Gross-Pitaevskii equation, the 
main result of this study so far. 

The hard-core coupling constant was found to be renor- 
malised by the presence of the random interactions in 
a temperature and spatially dependent way, though the 
spatial dependence of the renormalisation was seen to 
disappear in a mean-field limit where the variances of 
the Gaussian couplings are given by (t'^{x,x') = a'^/N 
and loose their spatial characteristics. The contribution 
from the random interactions to the equation for the con- 
densate contained the replica diagonal and replica off- 
diagonal Green's function for which the replica symmet- 
ric Dyson equations were given. These equations were 
formulated in terms normal and anomalous Green's func- 
tions. 

Having derived a set of equations to be satisfied by 
the condensate wave-function and the Green's functions, 
a next phase in this study is to evaluate these numeri- 
cally. Data for the inverse replica-diagonal Green's func- 
tion matrix may then be used to compute the dispersion 
relation for quasi-particle excitations. 

* * ^ 

This project has profited considerably from numerous 
discussions with R. Kiihn for which I am thankful. 
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